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The possibility of inducing topological superconductivity with cuprate high-temperature superconductors 
(HTSC) is studied for various heterostructures. We first consider a ballistic planar junction between a HTSC 
and a metallic ferromagnet. We assume that inversion symmetry breaking at the tunnel barrier gives rise to 
Rashba spin-orbit coupling in the barrier and allows equal-spin triplet superconductivity to exist in the ferro- 
magnet. Bogoliubov-de Gennes equations are obtained by explicitly modeling the barrier, and taking account 
of the transport anisotropy in the HTSC. By making use of the self-consistent boundary conditions and solu- 
tions for the barrier and HTSC regions, an effective equation of motion for the ferromagnet is obtained where 
Andreev scattering at the barrier is incorporated as a boundary condition for the ferromagnetic region. For a 
ferromagnet layer deposited on a (100) facet of the HTSC, triplet p-wave superconductivity is induced. For the 
layer deposited on a (110) facet, the induced gap does not have the p-wave orbital character, but has an even 
orbital symmetry and an odd dependence on energy. For the layer on the (001) facet, an exotic /-wave super- 
conductivity is induced. We also consider the induced triplet gap in a one-dimensional half-metallic nanowire 
deposited on a (001) facet of a HTSC. Due to the breaking of translational symmetry in the direction perpendic- 
ular to the wire axis, the expression for the gap receives contributions from different perpendicular momentum 
eigenstates in the superconductor. We find that for a wire axis along the a-axis, these different contributions 
constructively interfere and give rise to a robust triplet p-wave gap. For a wire oriented 45° away from the 
a-axis the different contributions destructively interfere and the induced triplet p-wave gap vanishes. For the 
appropriately oriented wire, the induced p-wave gap may give rise to Majorana fermions at the ends of the 
half-metallic wire. In light of the recent experimental progress in the quest for realizing Majorana fermions, we 
also discuss inducing superconductivity in the spin-orbit coupled nanowire using a HTSC. Based on our result, 
topological superconductivity in a semi-conductor nanowire may be possible given that it is oriented along the 
a-axis of the HTSC. 

PACS numbers: 74.45,+c, 74.72.-h, 74.78.Na 



I. INTRODUCTION 

The study of heterostructures involving a superconductor 
and non-superconducting materials has been a subject of in- 
terest for many years. However, the field has recently en- 
joyed a renewed level of intense activity since these sys- 
tems can harbor exotic emergent particles known as Majorana 
fermionsP^ These particles are sought after for their uncon- 
ventional non-Abelian anyonic statistics and for their poten- 
tial use in fault tolerant topological quantum computationP2J 
To this date, possible realizations of Majorana fermions have 
been discussed in a number of systems including 3 HeP frac- 
tional q uantu m Hall systemsp^ organic superconductors, 12 
SrRuC>3p^J fermionic cold atoms systems ESGH anc [ the sur- 
face state of a topological insulator! 1 ! 2 ! 18 ! 19 ! However, a series 
of theoretical works recently showed that Majorana fermions 
can be realized in superconductor heterostructures with rela- 
tively conventional building blocks ja eEQEl I The most promis- 
ing of these proposals suggests that a semiconductor nanowire 
fabricated on top of an s-wave superconductor, under the right 
conditions, should host a pair of localized Majorana excita- 
tions on the two ends of the wire! 4 ! 5 ! 20 -! The zero-bias con- 
ductance peak observed recently can be considered as an ev- 
idence for these Majorana excitations. 22 Very recently, trans- 
port properties compatible with the existence of topological 
superconductivity has also been observed in lead nanowires! 2 ^! 



The basis behind most of these recent Majorana propos- 
als is to "engineer" a single-component p x + ip y supercon- 
ducting state in a system which does not intrinsically exhibit 
this property. In th e semiconductor nanowire proposal men- 
tioned abovepEED a w jj- e w ith spin-orbit coupling is placed 
in contact with a conventional superconductor to produce an 
effective p x + ip v state in the wire. Upon applying a mag- 
netic field to the system, one of the components is effectively 
removed and the desired single-component p + ip topolog- 
ical superconductor is obtained. Another natural possibility 
is to induce superconductivity in a half-metallic ferromagnet. 
The system is metallic for one spin component while an insu- 
lator for the other and, as such, has a single non-degenerate 
Fermi surface. If triplet p-wave superconductivity is induced 
in such a system it must harbor Majorana fermions, either in 
vortices in a two-dimensional geometry! 11 ! 24 ^ at the ends in a 
one-dimensional wire geometry. 25 Indeed, proximity-induced 
topological superconductivity in a half-metallic ferromagnet 
in contact with a conventional superconductor has already 
been studied theoretically. 2 ^ 8 -! 

Superconducting proximity effect in metallic and half- 
metallic ferromagnets in contact with conventional su- 
perconductors has received much theoretical 2 ^ 5 -! an d 
experimental 4 ^" 5 - 5 -! attention over the years. It is known that 
in a uniform metallic ferromagnet singlet superconducing 
correlations penetrate over much shorter distances than in a 



non-ferromagnetic normal metal!^ 9 -! This is because singlet 
superconducting correlations compete with ferromagnetism 
which favors parallel alignment of electron spins. However, 
early transport data reported unexpectedly long-ranged 
proximity effect in Cc 46 and Ni 47 , where the decay length for 
superconducting correlations was typical of proximity effect 
observed in a normal metal. This long-ranged proximity 
effect was attributed to the formation of equal-spin triplet 
pairs which ar e imm une to the pair-breaking effects of an 
exchange field ! 30 ! 31 1 It was conjectured that local inhomo- 
geneity in the magnetization near the S/F interface induces 
triplet-pairing inside the ferromagnet, and generates an exotic 
odd-frequency s-wave symmetry for the triplet component 
of the condensate which should be robust to disorderPEEO 
For a half-metallic ferromagnet, singlet superconductivity is 
expected to be destroyed due to the vanishing density of states 
for the minority spin component. However, supercurrent flow 
was remarkably observed between two conventional super- 
conductors separated by CrC>2 (a half-metal) with a width of 
order 1/um. The supercurrent was attributed to odd-frequency 
triplet Cooper pairs which can be generated by spin-flip 
scattering at junction interfaces \2M2\ Singlet-to-triplet con- 
version mechanism due to both spin-mixing and spin-flip 
process, which is inevitably present when magnetic moment 
near the in terfac e is inhomogeneous, was also studied in 
this contextP^ISll Josephson coupling between a singlet and 
triplet superconductor in the presence of interface spin-orbit 
coupling was also theoretically studied. 57 In the context 
of topological superconductivity, theoretical works have 
discussed another scenario for singlet-to-triplet conversion 
based on spin-orbit coupling, which can exist e ither in the 
superconductoi^or at the S/F interface. ^ 28 * 58 * 59 ! 

The reliance on the superconducting proximity effect is 
paramount in almost all of the current proposals for realizing 
topological superconductivityP^inill] However, much of the 
recent interest has been in artificial creation of a simulated p- 
wave superconductors in heterostructures where the supercon- 
ductivity is proximity-induced by a conventional s-wave su- 
perconductor. In contrast to all of these past works, the current 
work considers the creation of topological superconductivity 
and Majorana fermions by using a cuprate high-T c supercon- 
ductor (HTSC) as the agent inducing the proximity effect. 
HTSCs are characterized by high transition temperatures, 60 
and may allow for more favorable experimental conditions 
for realizing Majorana fermions. In context with the pre- 
ceding discussion, we focus on the physics of proximity ef- 
fect between HTSCs and metallic (or half-metallic) ferromag- 
nets. Proximity effect in HTSC-ferromagnet and/or HTSC- 
half-metal heterostructures have been studied theoretically 
with a combination of extended Hubbard model, Hartree- 
Fock theory and Bogoliubov-de Gennes equationj^El m the 
presence of a spin-active interface, 64 65 with spin-conserved 
interface tu nnel ing and spin-bandwidth as ymme try in the 
ferromagnetpSEll and within circuit theorypS Enhanced 
proximity effect in the c-axis direction was also studied in a 
diffusive ferromagnet in the presence of a domain wall near 
the interface. 70 More recently, unconventional superconduc- 
tor/ferromagnet junctions on the surface of a topological insu- 



lator have also been investigated. On the experimental front, 
signatures of long-ranged Josephson coupling were observed 
for a YBa 2 Cu 3 CV,5 (YBCO)-SrRu0 3 (itinerant ferromagnet) 
junction. 72 More recent works also showed superconducting 
proximity effect in thin-film heterostructures comprising of 
the ha lf-me tallic manganite La2/3Cai/3MnC»3 (LCMO) and 
YBCO .! 48 ! 73 " 76 ! These experiments also observed long-ranged 
superconducting correlations in the hal f-meta l, implying spin- 
triplet pair formation in the manganite! 31 * 77 ! 78 ! 

Motivated by intense recent activity in the field of topo- 
logical superconductivity as well as significant experimen- 
tal progress mentioned above, we study proximity effect be- 
tween HTSCs and metallic ferromagnets. We first con- 
sider a ballistic planar junction composed of the two mate- 
rials separate d by a pot ential barrier containing Rashba spin- 
orbit coupling! 26 ! 28 ! 58 ! 59 ! The interfacial Rashba spin-orbit cou- 
pling provides the spin-flip mechanism with which equal-spin 
triplet pairing in the ferromagnet can be generated. We fo- 
cus on the proximity-induced mini-gap in the ferromagnet 
in the equal-spin triplet channel. In contrast to many of the 
past works, we do not employ the tunneling Hamiltonian for- 
malism to model the interface. Instead, the barrier region is 
resolved into a region of finite thickness a, and we micro- 
scopically model the interface with spin-orbit coupling. We 
study proximity effect in both the ab-plane and c-axis di- 
rections, taking account of the hopping interlayer transport 
in the superconductor. We consider three different crystallo- 
graphic orientations: (i) a-axis parallel to the interface normal 
((100) interface); (ii) a-axis 45° away from the interface nor- 
mal ((110) interface); and (iii) c-axis parallel to the interface 
normal ((001) interface). We also consider a ferromagnetic 
nanowire in proximity to a HTSC. Analytic expressions for 
the induced mini-gaps are provided for all the cases. The het- 
erostructure under consideration and the gap orientation cor- 
responding to each of the interface orientations are shown in 
Fig.0 

An important issue when inducing topological supercon- 
ductivity via proximity effect is the size of the induced mini- 
gap in the normal region. This mini-gap is crucial for the pro- 
tection of Majorana fermions from external perturbations. A 
motivation to use a HTSC in this context is due to the robust 
gap it offers even at temperatures of order 100K. However, it 
is possible that the induced mini-gap amplitude is limited by 
the energy scale set by the interface spin-orbit coupling, Asoc- 
This is indeed the case for orientation (i), where the mini-gap 
is independent of the gap amplitude for the proximate super- 
conductor. For this case, a larger gap in the superconductor 
therefore does not directly increase the size of the mini-gap. 
On the other hand, a proportional dependence on the gap am- 
plitude is found for the mini-gap in case (iii) (proximity effect 
in the c-axis direction). In this case, increasing the gap in the 
proximate superconductor can directly enhance the robustness 
of the gap in the normal region. We add that even when the 
mini-gap does not explicitly depend on the gap of the prox- 
imate superconductor, it is still essential for a robust gap to 
be present in the superconductor to have any kind of gap in 
the normal region. In this respect, both a large Asoc an d the 
proximity to a HTSC are a favorable combination for indue - 
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FIG. 1: (Color online) The tunnel junction considered. The ferro- 
magnetic metal, the spin-orbit coupled barrier and the superconduc- 
tor occupy < z < d, -a < z < and z < -a, respectively. Various 
crystal orientations for the HTSC are considered: (i) (100) parallel to 
the interface normal; (ii) (110) parallel to the interface normal; and 
(iii) (001) parallel to the interface normal. Gap with positive (nega- 
tive) sign is indicated by red (blue). The existence of an easy axis in 
the z-direction is assumed in the ferromagnet. 



ing a large mini-gap. Ref. [79] purports that certain metallic 
surface states can be subject to spin-orbit coupling strength 
of the order Asoc ~ lOOmeV. If spin-orbit coupling of this 
scale can be generated in the heterostructure considered here, 
having a HTSC should certainly be beneficial for generating 
a larger mini-gap in the normal region. An important exper- 
imental question in realizing Majorana fermions with HTSC 
is the scattering into the nodal quasiparticle states. This issue 
relates to the stability of the Majorana end states, and will be 
addressed in a forthcoming publication. 

The paper is organized as follows. The heterostructure 
and its theoretical modeling are introduced and discussed in 
Sec. [II] An effective Bogoliubov-de Gennes equation for the 
ferromagnetic region is obtained in Sec. [Ell] Analytical re- 
sults for the proximity-induced triplet mini-gap are discussed 
in Sec. IV for various crystallographic orientations. Induced 
triplet mini-gap for a half-metallic nanowire is studied in 
Sec. [V] A summary and experimental implications of our 



work are presented in Sec. VI 



II. THE SETUP AND MODEL 

We consider a ballistic planar tunnel junction between a fer- 
romagnetic metal (0 < z < d) and a HTSC (z < -a) as shown 
in Fig. [T] The two systems are separated by a tunnel bar- 
rier {-a < z < 0) containing 2D Rashba spin-orbit coupling, 
which can arise due to breaking of inversion symmetry about 
the interface plane. In the presence of spin-orbit coupling and 
superconductivity it is most natural to express the equations 
of motion in terms of a four-component vector in Nambu-spin 
space, [Vixj)] 7 = (iA T (*,r) t^CM) <ff\(x,t) <p\(x,t)), 
where x = (r,z) collects the coordinates r - (x,y) parallel 
to the interface plane and the transverse coordinate z- We note 



again that we do not employ the tunneling Hamiltonian for- 
malism to model the interface. Instead, the barrier region is 
resolved into a region of finite thickness a, and we microscop- 
ically model the interface with spin-orbit coupling. We there- 
fore start with equations of motion for three regions, each of 
which can generally be written as 



dx' [sf 3 6(x - x') - <H p (x,x')} y p {x',s) = 0, (1) 



where p - {F, B, S] labels the ferromagnet, barrier and su- 
perconductor regions and s is the Fourier variable for time, 
Wfat) = J ds/(2n) y ¥P(x,e)e- ist/h . We use f-, to denote the 
vector of Pauli matrices acting on Nambu space and "calli- 
graphic" letters, such as H, to denote 4x4 matrices acting on 
Nambu-spin space. Hats will also be used in general to denote 
2x2 matrices. 

The ferromagnet is described within the standard Stoner 
model, assuming identical dispersion for both spins and with 
a shift in energy by the exchange interaction h. Assuming that 
the ferromagnet has an easy axis along the z-axis, its Hamil- 
tonian is given by ( H F {x,x') = 5(x - x'YH F {x) with 



<H F {x) = f Q H F (x), 



(2) 



where H F {x) = -h 2 V 2 ./2mf — (Xf - mp and fif are effec- 
tive mass and chemical potential in the ferromagnet respec- 
tively, and <x, is used to denote the vector of Pauli matrices 
acting on spin space. 

The barrier in between the ferromagnet and the supercon- 
ductor is characterized by a large constant potential Uq and 
Rashba spin-orbit coupling with strength a. Its Hamiltonian 
reads e H B (x,x') = 6{x - x'YH B (x) with 



<H B (x) = 



H B (x) 
[H B (x)Y 



(3) 



where H B (x) = -% 2 Vl/2m B + Uo + a[& x (-iV x )] z , and m B is 
the effective mass in the region. 

Finally, we model the HTSC within the Bogoliubov-de 
Gennes mean-field theory and consider various orientations 
of its a-axis with respect to the interface normal. The Hamil- 
tonian can generally be written as 



9{ s (x,x') = 



6(x - x')Hl (x) A(jc - x') 
A*(x-x') 5(x - x')H?.(x) 



where 



A(x - x') 



A(x-Jt') 
-A(x' -x) 



(4) 



(5) 



The dispersion in the HTSC depends on its crystallographic 
orientation. We consider three different orientations where the 
interface is: (i) the (100) surface; (ii) the (110) surface; and 
(iii) the (001) surface (see Fig.[T|. For cases (i) and (ii), we 
use the bulk dispersion of the form 



J cos(-id x s) ■ 



2m s 



■Us, 



(6) 
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FIG. 2: (Color online) (a) 3D Fermi surface corresponding to dis- 
persion {6]l on the left and quasi-2D Fermi surface for the majority 
spin component in the ferromagnet on the right. An overlap between 
the two Fermi surfaces is achieved for arbitrarily small k^. (b) 3D 
Fermi surface corresponding to dispersion l|7j on the left and quasi- 
2D Fermi surface for the majority spin component in the ferromagnet 
on the right. An overlap between the two Fermi surfaces is achieved 
for JKL < JfcjL < JfcS. + 27. 



where J is the c-axis (interlayer) hopping amplitude, s is the 
distance between the copper oxide planes, and m$ and are 
the effective mass and chemical potential. For case (iii) (rel- 
evant for proximity effect in the c-axis direction) the relevant 
dispersion is 

#011 M = " ^ + J cos(-^i) - [i s ■ (7) 

The gap function in Q reflects the cZ-wave nature of the gap, 
and also depends on the crystal axis orientation. Its explicit 
expression will be given in the relevant sections below. 

We now discuss an assumption which will be made 
throughout the work (unless otherwise stated). It is an as- 
sumption regarding the magnitude of the Fermi wave-vector 
in the ferromagnet with respect to the Fermi wave-vector in 
the superconductor. The outcome of this discussion can po- 
tentially influence the experimental conditions favorable for a 
robust proximity effect between these two systems. For the 
clean planar tunnel junction considered in this work, momen- 
tum parallel to the interface must be conserved during tunnel- 
ing. This usually places a restriction on the region of the Fermi 
surface from which Cooper pairs can tunnelP^In Fig.j^a), the 
3D Fermi surface corresponding to dispersion (|6]l is plotted on 



the left and the 2D Fermi surface for the ferromagnetic metal 
is plotted on the right. This corresponds to the crystal orien- 
tations (i) and (ii) defined above (see Fig. [T}. In Fig. |2jb), 
the 3D Fermi surface corresponding to dispersion (|7]i is plot- 
ted on the left and the 2D Fermi surface for the ferromag- 
netic metal is plotted on the right. This corresponds to the 
crystal orientation (iii). We are considering a quasi-2D fer- 
romagnetic metal, which is confined in the z-direction, and 
focus our attention on the majority spin component since we 
are interested in inducing equal-spin triplet Cooper pairs. We 
assume that the Fermi wave-vector in the ferromagnet, k^, 
is controlled by externally gating the ferromagnet. Projec- 
tions of the Fermi surfaces (ring) of the ferromagnet onto the 
superconductor's Fermi surface are shown by the red lines. 
Cooper pairs which participate in the tunneling are expected 
to come from these regions because the in-plane momenta are 
conserved for these processes. For orientations (i) and (ii) (i.e. 
Fig.j^a)), the Fermi surface with an arbitrarily small k^ will 
have nonzero overlap with the Fermi surface of the supercon- 
ductor. In this case, a robust proximity effect is expected even 
when the ferromagnet is in the dilute limit (small k F F ^). On the 
other hand, for orientation (iii) (i.e. Fig. [2jb)), a finite over- 
lap with the Fermi surface of the superconductor is achieved 
only for ft™. < k^ < ft™ + 2J within our model. Therefore, 
ferromagnets with Fermi wave-vectors lying outside this re- 
gion may experience a strong suppression in the proximity ef- 
fect. This in principle places a restriction on the ferromagnet's 
wave-vector which is favorable for strong proximity effect. In 
the remainder of this work, we will assume a small ft™ <K k s „ 
for cases (i) and (ii), where H is the Fermi wave-vector in the 
superconductor, and ft™. < ftL < ft™. + 2J for case (iii). 

In the following sections, Andreev physics at a HTSC- 
ferromagnet interface in the presence of interface Rashba 
spin-orbit coupling will be studied. We note the reader that the 
formal procedure used to obtain the proximity-induced mini- 
gaps is applicable both to a polarized as well as unpolarized 
normal regions. However, the ferromagnetism in the normal 
region is crucial if it contains multiple transverse sub-bands 
and one wants to realize an isolated Majorana mode, i.e. (# of 
Majorana modes) mod 2 = 1. The usual difficulty in yielding 
an odd number of sub-bands at the Fermi level arises because 
the electrons have spin and consequently come in pairs. How- 
ever, an isolated Majorana mode is in principle obtainable if 
the normal region is a ferromagnet so that the total number of 
occupied sub-bands is odd, i.e. (A^ + AjJ mod 2 = 1, and if a 
p-wave equal-spin triplet superconductivity is induced in each 
sub -band. 



III. BOGOLIUBOV-DE GENNES EQUATION FOR THE 
FERROMAGNET 

The first goal is to obtain the effective equation of motion 
for the ferromagnet. This is done by first imposing proper 
boundary conditions at each of the two interfaces, then solv- 
ing the barrier and superconductor regions, and eventually in- 
corporating their effects as a boundary condition for the fer- 
romagnetic region. Due to translational symmetry in the xy- 



4 



plane (discrete translational symmetry in the x-direction when 
the copper oxide planes are stacked parallel to the yz plane) 
we may introduce a Fourier variable k = (k x , k y ) defined via 
^(z, e) = J d 2 r *P p (jc, e)e~' k r . Boundary conditions are then 
defined in the usual way. The first set of boundary conditions 
are constructed by imposing continuity of the wave-function 
at each interface, 

%(0, e) = y B k (0, s), y B k (-a, s) = W s k (-a, s), (8) 

and vanishing of the wave-function at the hard wall, 



%{d,s) 



0. 



(9) 



Boundary conditions for the derivatives are obtained by inte- 
grating the equations of motion over a small interval with re- 
spect to z and subsequently sending the interval to zero. Due 
to the particular two-dimensional nature of the spin-orbit cou- 
pling, the derivative boundary conditions are simply given by 



— d z w F k (o,s) = — d z v B k (o,s), 

nif Mb 



1 



1 

m s 



d^l(-a,s) = —d,y s k (-a,s). 



(10) 
(11) 



Recall that the equation of motion for the ferromagnet is 
given by ([T]i and Q. This equation of motion must be accom- 
panied by appropriate boundary conditions. At z — d, the hard 
wall boundary condition (|9]l applies. The boundary condition 
at z = is obtained by first solving ([T]) for the barrier and the 
superconductor and by making use of the boundary conditions 



<JsJ, ( 10 1 and ( 11 ipMJ] Incorporating this boun dary condition 
into the equation of motion for the ferromagnetpSEHwe obtain 



j dz'[g F k {z,z',e)Y l %{z',s) 



h 2 



m B 



Kk + T k 



H 2 /2m B 



1 - £- B er f (^- 



(12) 



where the inverse Green function [@ k (z,z',s)] 1 = ef^6(z - 
z')-'H F (z,z'),and<H F (x,x') = j d 2 kl(2n) 2 <H F k (z,z')e ik(r - r '\ 
We note that for cases (i) and (ii), continuous translational 
symmetry in one of the directions parallel to the interface 
plane is broken down to discrete symmetry associated with 
the periodic stacking of the copper oxide planes. In this case, 
the momentum Fourier variable corresponding to that spa- 
tial direction must be restricted to the first Brillouin zone, 
[-n/ s,n/s]. The transmission and reflection matrices read 



n 



t t o 

f 



-k 



R k 
R 



-k 



(13) 



Here, the asterisk denotes complex conjugation (without 
transpose) and the 2x2 transmission and reflection matrices 
are 



R k = 



-iK e 



(14) 
(15) 



We have also defined Kf - q+/(2 sinh(q + a))±q-/(2 sinh(g_a)) 
and k* = q+/(2tanh(q+a)) ± q-/(2 tanh(g_a)), with q ± = 

[(2mB/n 2 ){Uo-s + (h 2 k 2 /2m B )±ha\k\)] l/2 being the in- 
verse penetration length of a particle with a positive/negative 
chirality inside the barrier, and a is the barrier width. The an- 
gle jk is defined via k = k (cos y k , sin y^). Equation (12i is 
the Bogoliubov-de Gennes equation for the ferromagnet re- 
gion with superconducting correlations and the barrier spin- 
orbit coupling effects introduced via the boundary term. In 



( 14 1 and ( 15 1, the coefficient k\ is associated with tunneling 



of particles without spin flip and with tunneling of particles 
accompanied by a spin flip. Similarly, coefficients k* are as- 
sociated with the reflection of particles without (+) and with 
(-) a spin flip. As the width of the barrier grows, Kf decay 
exponentially, as expected. On the other hand, if the spin- 
orbit coupling in the barrier vanishes (i.e. a = 0), we see that 
the spin-mixing terms vanish, k~ - k~ - 0, and we recover 
boundary conditions for the standard proximity effectP 

In (jT2j, we have also defined the "surface" Nambu-Gor'kov 
Green function for the superconductor, which is evaluated 
at the interface between the barrier and the superconduc- 
tor, Q k mi (s) = \im z ^ a g s k (z,-a,e). Here, Q s k {z,z') is the 
Nambu-Gor'kov Green function matrix of the Bogoliubov- 
de Gennes equation for the superconductor valid on the half- 
space, z < -a. In the low-transparency limit, one may approx- 
imate Q s k {z, z') with a Green function for an isolated supercon- 
ductor. It then satisfies the Neumann boundary condition (zero 
flux condition at the right boundary), 



d:Q 3 Az,z')\ 



0. 



(16) 



The explicit form for the surface Green function matrix 
@ k mt (s) will depend on the crystallographic orientation. We 
will now compute the matrix next. 



A. Surface Green function matrix 

In this section, the surface Green function matrix @ k ud (s), 
which appears in ( p"2| ), is evaluated for the three different crys- 
tallographic orientations. The Green function matrix encodes 
information about the superconductor and crucially affects 
the Andreev physics which takes place at the junction. The 
general technique we employ to evaluate the superconductor 
Green function valid on the half^pace z < -a is the method 
of image and its generalizationP^IMl 



1. Proximity effect in the ( 100) direction 



We begin by obtaining the surface Green function matrix 
@ k mi (s) relevant for proximity effect in the (100) direction. 
Here, the copper oxide planes are stacked parallel to the yz 
plane and the a-axis of the superconductor is oriented parallel 
to the interface normal. In this case, the superconductor pos- 
sesses reflection symmetry with respect to the in terface plane, 
hence one can make use of the method of i ma g^S§|82|83| t0 ex _ 
press the Green function valid on the half-space in terms of 



5 



the infinite bulk Green function. Similar to previous worlP^ 
we have 



@l(z,z', e) = §l(z - z, e) + g\(z + z' + 2a, e), 



(17) 



where Q s k (z, s) is found by Fourier transforming the infinite 
bulk Nambu-Gor'kov Green function matrix 



Ql(z,e) 



Gk,k z ( s ) F kJk (s) 
F-k-kX-s) G-k,-kX-e) 



where 



Fk,kXs) = 
G k ,kXe) = 







Fk,kXs) 



-F- k ,-kX-s) 

G k ,kXs) 
G k ,kXs) 



(18) 

(19) 
(20) 



and the (time-ordered) normal and Gor'kov Green functions 
are defined as usual as 



Gk,k z ( E ) 
Fk,kX £ ) 



e+e 



k.k- 



e 2 -[^(*,y] 2 -A 2 A +*y 

e 2 -[^i(t,* z )] 2 -A2 



(21) 
(22) 



where rj «c 1 is a small parameter. Here, the dispersion is 
given by the form For this orientation the gap has the 
form 



k~ ky 



(23) 



We now invoke the assumption that was made in the discus- 
sion at the end of Sec. [II] Within this assumption, most of the 
electrons which tunnel into the ferromagnet are normal to the 
interface plane (i.e. \k x \ ~ \k y \ ~ k*p <K k s F ). Assuming also 
that we have strong anisotropy (i.e. J <k /z s ) we can then ap- 
proximate * h 2 k 2 /2ms —fis and A^. w Ao. Also assuming 
that |e| «: A , we find 



F k (z,s) 



I 



Aoe^ 



2n i^l \ 2 , • 



2 »v*" 



(24) 



Here, v^. is the Fermi velocity in the superconductor. Here, we 
will drop the normal component (which can be assumed small 
compared to the anomalous component). The local Green 
function matrix then is given by 



Sf 1 ^) * ^y?i0<5-2). 



(25) 



2. Proximity effect in the (110) direction 

We consider again the case where the copper oxide planes 
are stacked parallel to the yz plane, but this time rotate the a- 
axis of the superconductor by 45°. Here, the nodal direction 



points parallel to the interface normal. The gap then has the 
form 



k\+k 2 y 



(26) 



Since the superconductor breaks reflection symmetry with re- 
spect to the interface plane the method of image must be 
generalized r 2 l 83 IWe find that the following form for the Green 
function provides the solution, 

&l(z,z',s) =@ s k (z-z'-,s)-@ s k (z + z'- + 2a,s) 

x [d z g s k (a + z", s)]' 1 [d z g s k (-a - z'-, e)] , (27) 
- + . One can verify that the Neumann bound- 



16 1 is satisfied by (27 1. One can also see that 



where z?~ = z' 
ary condition 

the real Green function (the first term) satisfies the equation 
of motion for the Green function with a delta-function source 
located at z — zf~ where z'~ < -a. The regulator + > 
makes sure that the delta-function source is located in the 
appropriate half-space z < -a when z' = -a. The image 
Green function (the second term) on the other hand satisfies 
the equation of motion with the delta-function source appro- 
priately located outside the relevant half-space domain. One 
can also reproduce the method of image solution ( fT7| from 
( 27 1 for the case when reflection symmetry is present. In 
that case, the Nambu-Gor'kov Green function matrix obeys 
@\{z, s) = @ s k (-z, e). Then the derivatives must be odd, i.e. 
d z Q s Az,s) = -d z Q s A-z,s). Using this fact, En\ reduces to 



-a, the Green function matrix is then 



( 17 1 as expected. 

On the surface, z 
given by 

QT\s) = #£(0,e){i - [d z g s k (o-, s)]" 1 [d z g s k (o + ,e)]Y 

(28) 

With the same assumptions as for the (100) case, most of the 
electrons which tunnel through the barrier have small mo- 
menta in the in-plane direction. One may then approximate 
the gap (26 1 by A k ,k- ~ 2Aok z k y /(k s F ) 2 . We find for its normal 



component 



G k (0,e) 



I 



dkr 



I t, 2 k 2 \ 



2n 



(m x 



A A 9 



(29) 



where 5k = 2Aok y /kp, Due to the oddness of the k z integral, 
the anomalous component F^O, s) = 0. The derivatives of the 
Green functions are given by 



d z G k (0\s) 



I 



■i I , X-kl \ i0+z 
dk z lk z( e+ 2^-L>s)e 



2n 



f i h 2 k 2 \ a a 7 k 2 B 

^-(i^-^J -4A 2 ^+»7 



fop 
ftvf 



-d z G k (p-,e), 



(30) 
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and 



The surface Green function matrix then reads 



d z F k (0\s) 



S 



dk, 



ik z (-6 k k z /k s F )e i0+z 



In 



k' F 



+ ;?7 



S t F t (0-,c). 



Then inserting the above results, we obtain 



-2i 



F £ + 7T- + "7 



'6(e) 

&(6) <5 t 

<5* -&(e) 

-tf 4 



(31) 






-&(e)J 
(32) 



where % k {s) 



•('-£)■ 



3. Proximity effect in the ( 001 ) direction 

We now consider the case where the interface normal is par- 
allel to the c-axis. The gap function in this case depends only 
on the in-plane momenta and reads, 



At = A cos(2y t - 20 o ), 



(33) 



where 6q is the angle between the a-axis and the jt-axis. We 
will assume without loss of generality that 6q = 0. Due to the 
hopping nature in the z-direction, the dispersion which must 
be used here is given by (F7J, and the Fermi surface is given by 
a cylindrical shape as illustrated in Fig. [2jb). Here, reflection 
symmetry with respect to the interface plane is retained and 
the trivial method of image (c.f. (17i) can be used. In the 



strongly anisotropic limit J <c fig, the normal component for 
the surface Green function reads 



G k (z = Q,e) 



Tils 



' nls dk 7 



" 2 " e 2 • 



A 2 k + it] 



s e 2 



A\ + iTJ 



(34) 



where the in-plane dispersion = fr 2 k 2 /2m s - (1$. For the 
anomalous component, 



F k (z 



X* dk 



1 + 



(35) 



QT'is) 



2/s 



r*+3 




Aa- 



A^ + 1?7 


-At 






A* 




-A ft 





(36) 



The appropriate surface Green function matrix for a given ori- 
entation is inserted into (12i. By evaluating the right-hand 



side of the equation, we extract the induced triplet gap in the 
ferromagnet. The results are presented and discussed in the 
next section. 



IV. PROXIMITY-INDUCED TRIPLET GAP 

We now obtain analytical expressions for proximity- 
induced gaps in the ferromagnet, and focus on the equal- 
spin triplet pairing channel. In order to evaluate the effects 
of the Andreev processes at the interface, we must evaluate 



the right-hand side of ( 12 1. The details of the superconductor 
are now fully encoded in the surface Green function matrix, 
(e), which are evaluated for different crystal orientations 



in Sec. Ill A We begin by writing the ferromagnet state vector 



in a separable form 



^(-e) y^z), (37) 



where ifi (tl '\z) is the transverse envelope wave -function. We 
have ignored the possible spin-dependence of this transverse 
wave-function since we assume h to be much smaller than the 
Fermi energy scale in the ferromagnet, i.e. h <K fX F . A dis- 
cussion on obtaining (A (tl) (z) is provided in AppendixjA] Here, 
we consider a thin enough quantum well such that most of 
the particles occupy the lowest transverse sub-band character- 
ized by wave-function \jA f (z) and energy e'J. ~ ft 2 7r 2 /2m F d 2 



We then introduce solution ( 37 1 into ( 12 1, multiply it from the 
left by ifr^\z), and integrate it over the thickness of the ferro- 
magnet. This leads to an effective equation of motion for the 
majority spin component of the ferromagnet which generally 
has the form 



E g (s, k)(f>] 



(-£) = 0. 



(38) 



Here = H 2 k 2 /2m F - /j f is the 2D dispersion for the ferro- 
magnet, E g (s, k) is the proximity-induced triplet gap, which in 
general can have a frequency dependence, and we dropped a 
weakly induced spin-orbit coupling term and the singlet pair- 
ing term. 

The task now is to obtain this proximity-induced triplet gap 
Eg(s, k) for the various crystal orientations of the supercon- 
ductor. Crystal axis orientation has been known to play an 
important role, for instance, in the tunneling spectra for a nor- 
mal metal -HTSC heterostructureP^K^The proximity-induced 
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triplet gap in our context also shows a strong orientation de- 
pendence. The expression for the gap is obtained from the 



right-hand side of ( 12 1. In computing the matrix, we keep the 



spin-orbit coupling strength a finite only in the spin-mixing 
terms'^ The reflection and transmission matrices have been 
defined in ( fT3| l and the surface Green functions for various 
crystal orientations are evaluated in Sec. Ill A For proximity 
effect with the (100) surface (case (i)), the gap function has 
d x i-f symmetry and is given by A^ 1 " ' = Ao{k 2 - k 2 )/(k 2 + 
The induced triplet gap is given by 



(39) 



where / = 1/ftv*, (3 = m S K + /m B k s F and k = [^ r) (Q)f. 
Here, v s F and k s F denote Fermi velocity and Fermi wave- 
vector for the superconductor, respectively. For weak spin- 
orbit coupling (i.e. tiak F ^ <K Uq), k + » q/ tanh(ga) where 

q ~ ^2m B Uo/fi 2 . Thus, for a high enough barrier, i.e. qa > 1, 
k + w q. We will assume here that the barrier is large but its 
scale is still smaller than the Fermi energy scale in the su- 
perconductor, i.e. q <sc ki. In this case, j3 <K 1 (as long as 



m B ~ m s ) and ( 39 1 can be approximated as 



7(100) 



(*) 



He-^fkoh 4 



(40) 



We see from Eq. (40 1 that the proximity-induced gap has the 



p-wave orbital character which arises due to the presence of 
spin-orbit coupling in the barrier. Comparing this result with 



the results of Ref. 28 ( 40 1 is essentially identical to the triplet 
gap induced by an isotropic s-wave superconductor. This is 
because for k F ^ <k k s F (see discussion at the end of Sec. |5j , 
most of the electrons which tunnel have momentum directions 
that are nearly perpendicular to the interface plane. For the 
interface with the (100) surface these electrons are moving 
in the anti-nodal direction and are thus subject to the maximal 
value of the gap. In this sense, it is reasonable that the induced 
triplet gap here is very similar to the proximity-induced gap 
with an s-wave superconductor. 

We also note here that the magnitude of the induced gap 
is sensitive to the actual width of the quantum well. This 
sensitivity comes from its dependence on the value of the lo- 
cal transverse wave-function ko, which is given in (A3 I. We 



find that for smaller widths Andreev reflections occur multiple 
times and the gap is consequently enhanced. 

We now move on to the case where the interface is made 
with the (110) surface. The gap function here has d XY symme- 
try and is given by A^ 0) = 2Aok z k y /(k% +ky). This situation is 
very different from the previous case. Here, most of the tun- 
neling electrons have momenta in the nodal direction of the 
gap. We, therefore, find a suppression in the triplet gap in- 
duced on the ferromagnetic side. The gap for this case reads 
(again for ft <k 1) 



7(110) 



(«,*)*- 



2e-<^fkoh 4 



m-j. 



(41) 



where <5^ v = 2Aok v /k s F . We see that for electrons tunneling 
near normal to the interface (i.e. k y /k s F <K 1) the gap magni- 
tude is strongly suppressed. Furthermore, the gap symmetry 
here is not of the p-wave form. Apart from the p-wave factor 
e~' 7k , there is an additional factor 8^ , This factor is spatially 
odd with respect to the y-direction. This is because for the 
(110) interface the gap function is odd with respect to reflec- 
tions in the y-direction. This is intimately tied to the existence 
of zero-bias conductance peaks mentioned earlier. This ad- 
ditional factor makes the gap spatially even, but the required 
oddness is retained by the linear frequency factor making this 
gap odd in frequency. We have also computed the singlet gap 
for this interface orientation induced in a (non-ferromagnetic) 
normal metal and for zero spin-orbit coupling. We find that 
the induced singlet gap is also odd in frequency. However, 
an even-frequency gap results for the (100) surface interface. 
There may be a curious connection in this odd-frequency ex- 
pression for the gap and odd-frequency pairing which has 
been discussed extensively in t he context o f inhomogeneous 
superconductor heterostructuresPsEHEHH 

We now move on to proximity effect in the c-axis direction 
where the interface is parallel to the (001) surface. The gap 
function is given by At = Ao cos 2jk, and we may assume 
without loss of generality that the anti-nodal direction is along 
the jc-axis (see Fig. [TJ. Here, due to the hopping transport in 
the z-direction we do not assume k^ <sc Id,. As discussed at 

it is most natural to assume that k F ^ ~ k s F ' 2D 
ror maximal proximity effect. For the induced triplet gap, we 
find 



4 001) (e,it) 



(S) + MS)^-^-^) 2 - A *]' 



where £f = h 2 k 2 /2ms — fis ■ In the large barrier limit, 



qs 



k, £ 



U 



«: 1, 



142b can be approximated by 



7(001) 



(£,*) 



2ie- iY "sk h k K ^- 
q 2 



(42) 



(43) 



(44) 



We can then give an estimate for the product qs in ( |43] l. The 
interlayer distance between copper oxide planes can be esti- 
mated as s ~ A\ciq, where A\ is a number of order 1 - 10, 
and ao is the square lattice constant for the copper oxide 
planes. But since we are assuming that the barrier energy is 
smaller than the Fermi energy scale in the superconductor, i.e. 
q ~ k s F /A2 with Az > 1, the product qs can be estimated as 
qs ~ AxQdlA^a^ = A\jAi ~ 0(1). Therefore, the large bar- 



rier limit (43 i can be simplified to A k , s) <K Uq <s fis- 



From (44 1, we see that the induced triplet gap does not have a 
p-wave form but has an /-wave form proportional to Ajte~' r * . 



We note here that the induced mini-gap (44i is directly 
proportional to the gap amplitude Ao in the proximate super- 
conductor. This contrasts with the (100) case where Ao does 
not explicitly enter the final expression for the mini-gap (c.f. 
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Half-metallic nanowire 




High-T c Superconductor 
(b) 




FIG. 3: (Color online) (a) Half-metallic nanowire placed on the 
(001) surface of a HTSC. The copper oxide planes are indicated in 
light yellow, (b) The nanowire placed on the superconductor with 
different orientations with respect to the underlying gap symmetry. 
The x-axis is always defined along the axis of the wire. 



(40 1). Therefore, for the (001) case, a larger gap in the proxi- 
mate superconductor can directly enhance the mini-gap in the 
normal region. 



V. TRIPLET SUPERCONDUCTIVITY IN A 
HALF-METALLIC NANOWIRE 

We now consider placing a half-metallic nanowire on the 
(001) surface of a HTSC and discuss the equal-spin triplet 
mini-gap induced in the wire. As shown in Fig. [3ja) and 
(b), we envisage placing the nanowire in different orientations 
with respect to the underlying li-wave gap orientation. We de- 
fine the angle between the wire axis and the a-axis with 6q, 
and consider 9o = 0, n/4, n/2. As shown in Fig. [3] we will 
always take the jc-axis to be the direction along the axis of the 
wire. Due to the tunneling transport in the c-axis direction the 
physics of proximity effect here should not depend strongly 
on the number of copper oxide planes considered. Therefore, 
we simplify the problem by ignoring the coupling of the first 
copper oxide layer to the rest of the layers below, and consider 
the problem of a half-metallic nanowire deposited on a 2D su- 
perconductor with a af-wave order parameter. In contrast to 
the discussion at the end of Sec. [II] we assume that the wire is 
in a dilute limit where <sc k s F , where k s F (k^) is the Fermi 
wave-vector for the superconductor (majority spin species in 
the wire). 

We study the equal-spin triplet gap that is induced in 
the half-metallic nanowire within the path-integral approach. 
Since the only focus is on obtaining an analytic expression for 



the triplet mini-gap the half-metallicity of the nanowire is not 
strictly germane, so we relax the half-metallic condition and 
consider a ferromagnetic nanowire. The Matsubara action for 
the ID nanowire and the 2D superconductor can be written in 
the Nambu-spin basis as 



where the 4x4 inverse Green function matrix for the super- 
conductor is given by 



-ie n + hf k + (^2)('0"2)A t , 



(46) 



and e„ = (2n + V^nksT is the usual fermionic Matsubara 
frequency. The Green function matrix for the ferromagnetic 
nanowire is given by Q^{e„), but we refrain from writing 
down an expression for it explicitly since it will not be impor- 
tant in the rest of the discussion. For Oq = 0, the gap function 
is given by A° = A (k 2 x - k 2 )/(k 2 x + fy, while for O = n/4, it 

reads A* /4 = 2A k x k y /(k 2 x + k 2 ). We now define the tunneling 
Hamiltonian. The tunnel matrix elements will be constructed 
based on the results from Sec. [TTT] However, unlike Sec. [Ill] 
we assume that the barrier Rashba spin-orbit coupling resides 
only along the ID boundary region between the wire and the 
superconductor. Nevertheless, the presence of spin-orbit cou- 
pling in the barrier will give rise to both spin-flip terms and a 
momentum-dependence in the tunneling amplitudes. The tun- 
neling matrix elements should reflect the same symmetry as 
the ID Rashba spin-orbit coupling present at the boundary.^ 
We thus take them to be proportional to the transmission ma- 



trix given by ( 13 1 and ( 14 1 



-rlD _ 



f TT it u sgn(£ v ) 
-it u sgn(^) t u 



(47) 



and the tunneling Hamiltonian can be written as 



>[ o _pn- K(e„) + c.r.. 



(48) 



The total action is then given by S — So+S j. At this point, we 
integrate out the superconductor fields and obtain an effective 
action for the wire, 



(2;r) 



(49) 



The k y -integral over the superconductor Green function ma- 
trix, Q s k {e„), can now be computed. We will now focus on 
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FIG. 4: (Color online) The 2D Fermi surface of the superconduc- 
tor (thick dashed line) and the rf-wave gap are schematically plotted. 
Red (blue) region corresponds to the positive (negative) part of the 
gap. The Fermi wave-vector in the superconductor is given by k s F . 
For a dilute wire (i.e. k^ <k fc£), the momentum of a typical particle 
in the wire is also plotted for 9q = and 6q = n/4 (the short arrows). 
The x-axis is always defined to be along the axis of the wire. Due 
to translational symmetry breaking in the y-direction, momentum in 
that direction is not conserved during tunneling. In contrast, due 
to translational symmetry in the x-direction, k x must be conserved 
during tunneling. The dashed (dotted) line represents regions in mo- 
mentum space which can be explored by a particle as it tunnels from 
the wire into the superconductor. For 8 = regions in momentum 
space which can contribute to the induced gap in the wire is colored 
in green. For 8 = n/4, it is colored in yellow. 



the matrix component in the second term of (49 1 which gives 
the equal-spin triplet mini-gap. For 9q = 0, the triplet pairing 
component reads 



E»(k x )*2it n t u fsgn(k x ), 



(50) 



where / = \/tiv s F . For 6q = n/4, the integrand for the triplet 
pairing component becomes odd in k x and so we have E n J A - 
0. 



We see from ( 50 1 that the orbital symmetry of the pairs 



formed in the ferromagnet is determined by the tunneling am- 



plitudes defined in (47 1. The pairing occurs within this model 
because electrons with momenta k x and —k x in the wire hop 
into the superconductor, where they form a pair. Because 
of translational symmetry in the x-direction the momentum 
component k x must be conserved during hopping. On the 
other hand, due to translational symmetry breaking in the y- 
direction a pair of electrons with momenta k x and —k x can 
form a pair inside the superconductor with many k y values, 
i.e. Cooper pairs with momenta (k x , k y ) and (—k x , —ky) where 
ky can take on many values. The effective gap produced in the 
ferromagnet is found by summing over all of these possible 
Cooper pairs. 

In Fig. [4] the Fermi surface of the superconductor and the 
cZ-wave gap are schematically plotted. Assuming that we are 
at low temperatures, T <K Ao, an electron tunneling into the 
superconductor from the wire has momentum k x ~ k^. For 
6*o = 0, the orientation and amplitude of k x with respect to the 
Fermi surface of the superconductor is also plotted in Fig. [4] 



and labeled as k x . For 8q = n/4, it is labeled as k x 14 . The 
dashed and dotted lines are used to represent momenta a par- 
ticle can explore as it tunnels into the superconductor from 
the half-metal. The dashed line is for 6q — and the dot- 
ted line is for 0q = n/4. For 0q = 0, the summation over k y 
involves regions in momentum space where the gap has the 
same sign (green line segments in Fig. [4]). In this case, the dif- 
ferent regions give constructive interference and a robust gap 
is obtained. On the other hand, for 9q = n/4, the summation 
over ky involves regions in momentum space where the gap 
has equal but opposite sign (yellow line segments in Fig. ffl. 
In this case, the two regions destructively interfere and the gap 
vanishes. With the same token, strong proximity effect is ex- 
pected again for 9q = n/2, but the induced gap will have a 
sign opposite to the 9q — case. For 9q — and O = n/2, 
the induced gap in the wire has the p-wave orbital character 



(see ( 50 1) rendering the realization of Majorana fermions on 
the ends of the half-metallic wire possible!^ 

Now let us consider the case where the half-metallic wire 
is replaced by a non-ferromagnetic nanowire but with spin- 
orbit coupling. We now remove the spin-orbit coupling in the 



barrier (i.e. = in (47 1), and the resulting induced gap in 
the wire should have an s-wave orbital character. For a wire 
with the 9q = orientation the induced s-wave gap has the 
form 



4 inglet fe) = -2% 



hif> 



(51) 



while the gap vanishes again for 9q = n/4. We note that the 
9q = wire has the essential ingredients for realizing Majo- 
rana fermions on its ends. 



VI. SUMMARY AND EXPERIMENTAL IMPLICATIONS 

In this work, we study proximity effect between a HTSC 
and a metallic ferromagnet. We first consider a quasi-2D fer- 
romagnet deposited on (100), (110), and (001) facets of the 
superconductor. The induced equal-spin triplet mini-gap for 
each of the geometries is analytically evaluated. We assume 
that the breaking of mirror symmetry at the boundary region 
between the superconductor and ferromagnet gives rise to 
Rashba spin-orbit coupling in the region. In contrast to many 
of the past works, we do not employ the tunneling Hamilto- 
nian formalism to model the interface. Instead, the barrier 
region is resolved into a region of finite thickness, and the in- 
terface with spin-orbit coupling is microscopically modeled. 
We also consider a half-metallic nanowire placed on top of 
a (001) facet of a HTSC. The induced triplet mini-gap in the 
wire is obtained for various wire orientations with respect to 
the underlying orientation of the d-wave order parameter. We 
now highlight some relevant transport and proximity effect ex- 
periments involving high-T c cuprates, and then comment on 
some experimental implications based on the findings in this 
work. 

The in- and out-of-plane transport anisotropy is also known 
to affect the proximity effect involving HTSCs. Experimen- 
tally, proximity effect was found to be weaker in the c-axis di- 
rection than in the ab-plane direction due to a much shorter co- 



10 



herence length in the former direction! " 9 ' We note, however, 
that so-called giant proximity effect has been reported both 
in the ab-plana^H2ll and c-axis^ directions. Proximity effect 
in HTSC-ferromagnet and/or HTSC-half-metal heterostruc- 
tures have also been studied experimentally. Early studies 
showed long-ranged Josephson coupling for a YBCO-SrRu03 
(itinerant ferromagnet) junction, 72 though the results may be 
controversial. 97 More recent works showed superconducting 
proximity effect in thin-film heterostructures comprising of 
the ha lf-me tallic manganite La2/3Cai/3Mn03 (LCMO) and 
YBCOP^HEH These experiments also observed long-ranged 
superconducting correlations in the half- metal, imp lying spin- 
triplet pair formation in the manganite,' 30 ! 31 ! 77 ! 78 ! b u t further 
verification may still be neededpS A very recent work has 
shown evidence for MacMillan-Rowell resonance in the same 
LCMO-YBCO heterostructure, further substantiating the evi- 
dence for long-ranged proximity effect in this systemP^ 

We first comment on heterostructures built from a 2D 
metallic ferromagnet (or half-metal) deposited on one of the 
facets of a HTSC. Here, our results are likely to be applica- 
ble when the Fermi wave-vector for the majority spin species 
in the ferromagnet is much less than the Fermi wave-vector 
of the superconductor (see discussion at the end of Sec. [EJ. 
We find that a robust spin-triplet p-wave superconductivity is 
induced in a 2D half-metal if the layer is deposited on the 
(100) facet of the HTSC (c.f. gOj). On the other hand, de- 
positing the same layer on the (110) facet is expected to be 
less favorable from the point of view of obtaining a robust gap 
(c.f. (|4l}). This is because most of the electrons which tunnel 
from the superconductor into the ferromagnet have momenta 
in the nodal direction of the gap. As such, these electrons are 
not subject to the gap in the superconductor and a weaker gap 
is therefore induced in the ferromagnet. Furthermore, in this 
orientation, the induced gap does not have a p-wave orbital 
symmetry, but has an even orbital symmetry. The required 
overall oddness is maintained via an odd-dependence on en- 
ergy. For a half-metal layer deposited on a (001) facet, the sit- 
uation differs from the other two cases. In this case, a robust 
proximity effect is expected when the Fermi wave-vector in 
the half-metal is approximately equal to the (2D) Fermi wave- 
vector for electrons in a copper oxide plane (see discussion at 
the end of Sec. Iffl. The induced gap here has an /-wave or- 
bital symmetry. Therefore, from the point of view of obtaining 
a robust p-wave mini-gap, depositing the ferromagnetic layer 
on a (100) facet of the superconductor is the most favorable 
option. 

The expressions for the mini-gaps obtained in Sec.|IV|have 
been assumed to be uniform in the normal region due to strong 
transverse confinement. For a clean, thin normal region, a ro- 
bust s-wave singlet mini-gap co-exists with the triplet gap for 
the (100) interface, and is given by E™ s = (2fk n 4 /m 2 B )(4) 2 
(we obtain no S 2 = triplet mini -gap within our model). 
These mini-gap expressions are essentially valid also for an 
unpolarized metal since the existence of the triplet gap relies 
only on the presence of the Rashba spin-orbit coupling at the 
interface. However, as mentioned at the end of Sec. [TTJ the 
ferromagnetism in the normal region becomes crucial for re- 
alizing an isolated Majorana zero mode at the boundary of the 



region. 

One may also consider depositing a half-metallic nanowire 
on top of a HTSC. From the experimental point of view, per- 
haps the most natural geometry is to place the wire on top 
of a (001) facet of the superconductor. This is indeed the 
orientation we consider here, and we show that inducing an 
equal-spin triplet p-wave gap in the wire is indeed possible. 
However, one must note that the p-wave gap will not be in- 
duced in the wire when it is oriented along the nodal direction 
of the underlying gap, i.e. the wire axis is 45° away from the 
a-axis (see Fig. [3] and discussion in Sec.|V]l. The absence of 
the gap is due to destructive interference among different hop- 
ping processes which contribute to the proximity effect. For 
the wire orientation along the a- or the fe-axis, the different 
hopping processes constructively interfere and a robust gap is 
expected. 

We also add that this setup used for the half-metallic 
nanowire can also be applied to the case where a non- 
ferromagnetic spin-orbit coupled nanowire is placed on top 
of a (001) facet of the HTSC. This setup is more akin to the 
recent semi-c onduc tor nanowire proposal for realizing Majo- 
rana fermions! 4 ! 5 ! 22 ! We may apply the above nanowire study 
to the case where the Rashba spin-orbit coupling is introduced 
in the nanowire instead of the barrier region. In this case, for 
a wire oriented in the a-axis direction, an s-wave gap should 
be induced in the wire (see discussion in Sec. W). The gap, 
however, should vanish (or be strongly suppressed) when the 
wire is oriented 45° away from the a-axis. The system consid- 
ered here can serve as a starting point for the recent Majorana 
proposal presented in Refs. |4] and |5] An important experi- 
mental question in realizing Majorana fermions with HTSC 
is the scattering into the nodal quasiparticle states. This issue 
relates to the stability of the Majorana end states, and will be 
addressed in a forthcoming publication. 

While our discussion for the nanowire geometry was re- 
stricted to a strictly one-dimensional wire, recent theoretical 
works have shown that this one-dimensionality condition is 
not required for realizing Majorana fermionsP^KS] j n R e f s 
l99land flOOl it was shown that Majorana end states can still 
be realized if the wire width does not greatly exceed the su- 
perconducting coherence length and an odd number of trans- 
verse sub-bands are occupied. Furthermore, Refs. 1 1 1 1 and 
I102l studied how the stability of the topological superconduct- 
ing state can be enhanced in multi-band wires. 

Our findings suggest that for suitably designed hybrid struc- 
tures and with good contact established at the interfaces Ma- 
jorana fermions can be realized at the ends of either a ferro- 
magnetic or a spin-orbit coupled semi-conductor nanowire in 
contact with a HTSC. If a vortex can be created in a 2D geom- 
etry, a Majorana fermion will reside at the vortex core while 
the other will reside at the outer edge. The manipulation of 
this state is less well-studied but still interesting. 
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The boundary value of the transverse wave function can be 
obtained from jAl\ and reads 



Appendix A: Obtaining the transverse wave-function 

In this appendix, we evaluate the transverse wave-function 
for the particles residing in the quasi-2D ferromagnet. This 



quantity is used in ( 37 i and appears in all the expressions for 
the induced mini-gap. We will assume here that the exchange 
energy scale is smaller than the Fermi energy, h <sc fi F . In 
principle, the transverse wave-functions for the majority and 
minority species are different. For instance, in the limit where 
/) » //f, we expect that the transverse wave-function for the 
minority component will be an evanescent wave which decays 
with the scale {/, * [2/«p h/h 2 ]^ 1 ^ 2 . For h <k fif, however, the 
dependence can be ignored. 

The transverse wave-function can be approximated as a so- 
lution to the free ID Schodinger equation, 



2m f 

with the boundary conditions, 



h 2 d 2 

— ^<A (tr) (z) = s n ^ tr \z), 



^cL<A (tr) (0) = *V tr) (0), ^ lx) {d) = 0. 
nip 



(Al) 



(A2) 



k = [< r) (0)f 



1 



(A3) 



where £ = m / r^/(m B fe tr ) and k tr is a solution to the eigenvalue 
problem, tan (k a d) / (k a d) = -ms/inifgd), which determines 
the spectrum, e tl - = h 2 k 2 /(2mf). For the lowest sub-band so- 
lution, iff?' (z), the energy is given by ejj. = n 2 ^) 2 /(2mp), 
where k®. is the wave-vector corresponding to the lowest en- 
ergy sub-band. 



We note that k^ is proportional to 1/d where d is the width 
of the transverse confinement. As the width decreases ko in- 
creases, signifying the fact that the particles reflect back to the 
interface more frequently and, thus, has a higher probability 
of residing there. This amplifies the Andreev physics at the 
boundary and leads to the enhancement of the induced mini- 
gap. This was briefly discussed in Sec.|IV| 
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